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1.  INTRODUCTION 


Most,  if  not  all,  measurements  of  a physical  quantity,  involve  the  comparison 
of  a reference  quantity,  of  known  magnitude,  with  the  unknown.  As  a result  of 
the  comparison,  the  reference  quantity  is  adjusted  in  magnitude,  and  the  process 
repeated  until  the  comparison  shows  negligible  error.  The  magnitude  of  the 
reference  quantity  is  then  the  estimated  value  of  the  unknown. 

If  the  unknown  is  varying  with  time,  it  may  be  that  the  variation  within  one 
intercomparison  period  is  sufficient  to  prevent  the  achievement  of  a satisfactorily 
small  comparison  error  with  simple  adjustment  of  the  reference  quantity;  that  is 
the  system  suffers  from  a lag  error.  In  this  case,  an  improved  estimate  may  be 
obtained  if  the  rate  of  change  of  the  unknown  is  also  estimated,  and  an  additional 
adjustment  is  made  to  the  value  of  the  reference,  equal  to  the  product  of  the 
estimated  rate  of  change,  and  the  intercomparison  period. 

This  process  can  be  extended  to  include  an  estimate  of  higher  order  derivatives, 
each  one  being  used  to  adjust  the  estimate  of  the  lower  ones.  Apart  from 
increased  complexity,  penalties  include  the  greater  sensitivity  to  random  errors 
occurring  in  the  comparison  process,  and  the  need  for  a succession  of  comparisons 
before  reasonable  estimates  are  obtained  of  the  higher  derivatives. 

The  process  described  above  is  often  called  a sampled  data  tracking  system. 

It  breaks  down  logically  into  two  parts:  the  discriminator,  which  compares  the 
last  estimate  with  the  variable,  and  estimates  the  correction  required;  and 
the  tracking  filter^predictor , which  remembers  the  last  value  of  the  estimates, 
and  calculates  the  improved  estimates  for  the  next  comparison.  The  whole 
comprises  a sampled  data  feedback  controller,  which  exhibits  the  usual  stability 
problems,  and  transient  responses  to  initial  conditions  and  changes  in  the 
primary  variable. 

In  this  report,  the  third  order  exponential  tracking  filter-predictor,  and 
the  system  incorporating  it,  are  defined  and  analysed  in  some  detail,  to  determine  | 
the  basic  characteristics  of  stability,  transfer  function,  transient  response,  j 
tracking  lags,  and  sensitivity  to  discriminator  noise.  Lower  order  systems  are 
also  examined.  A brief  analysis  along  the  same  lines,  is  made  of  a third  order, 
continuous  data,  tracking  filter-predictor,  an  example  of  which  is  the  electro- 
mechanical range  or  bearing  tracker,  commonly  used  in  older  radar  systems.  It 
is  shown  that  the  characteristics  of  the  sampled  and  continuous  data  systems 
are  almost  identical,  providing  the  inter-sampling  period  is  much  less  than  the 
effective  time  constants  of  the  tracking  filter.  An  equivalence  table  of 
system  parameters  is  given. 

Familiarity  with  matrix  notation,  and  with  classical  control  theory  is  assumed. 
Statistical  parameters  are  derived  from  first  principles. 

2.  MODEL  DEFINITIONS  AND  PREDICTION  MATRIX 

A sampled  data  tracking  loop,  using  a third  order  filter-predictor  is  illus- 
trated in  figure  1.  The  variable,  x,  is  the  unknown,  and  x,  x and  5 are  the 
estimates  of  x and  its  first  two  derlvaties  with  respect  to  time.  The  suffixes, 
t-1,  t,  t+1,  denote  successive  sampling  (comparison)  times,  separated  by  the 
inter- comparison  period,  T.  The  double  suffix,  t/t-1,  denotes  an  estimate  for 
the  value  at  time,  t,  calculated  on  the  basis  of  data  available  at  time,  t-1, 
i.e.  a prediction. 

At  time,  t,  there  exists  an  (unobservable)  correction,  x^^,  which,  if  added 

to  the  predicted  value,  x^,  will  give  the  true  value,  x^.  The  correction  is 

mixed  with  noise,  or  other  random  errors  in  the  Physical  System  under  examina- 
tion, which  may  also  add  non-linearities  and  dispersions.  The  Correction 
Estimator  produces  an  estimate,  x , of  x which  contains  a noise  ccxnponent. 
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8^,  and  a multiplicative  error  coefficient,  D,  which  is  a fimction  of  the  correct- 
ion required.  This  coefficient,  D,  accounts  for  non-linearities  and  gain  errors 
in  the  overall  Discriminator  characteristic.  It  ideally  takes  the  value,  unity, 
but  always  exhibits  an  S-type  characteristic,  due  to  the  finite  range  of  sensitivity 
of  the  Discriminator.  (This  finite  range  is  necessary  if  two  or  more  physical 
quantities  must  be  resolved,  e.g.  two  radar  echoes,  closely  spaced.) 

A Sampling  Device  is  located  somewhere  within  the  Discriminator  section,  and 
is  shown  ahead  of  the  Correction  Estimator  in  figure  1.  In  the  case  of  a pulsed 
radar,  this  function  is  provided  by  the  pulse  itself.  Its  purpose  is  to 
staticise  the  operation  of  the  Estimator,  at  values  appropriate  to  the  sampling 
time,  ignoring  changes  in  variables  between  sampling  times. 

The  Tracking  Filter  combines  the  estimates  of  the  physical  quantity  and  its 

derivatives,  x,  x,  applicable  at  time,  t,  but  predicted  from  data  from  time, 
t-1,  with  the  estimated  correction,  x^^,  applicable  at  time,  t,  to  produce  a 
revised  estimate  of  x,  x and  i at  time,  t,  and  based  on  the  extra  data  available 
at  that  time  (i.e. 

The  Tracking  Predictor  carries  out  the  simple  extrapolation  to  predict  the 
values  of  x,  x,  and  X,  at  the  next  sampling  instant,  T seconds  away,  and  based 


on  information  at  time,  t. 
Thus : 

\*l/t 

= T.  x^^^  + . x^^^ 

> 

(1) 

Vi/t 

= 0.5  T?  ^x^/^  ^ 

Note  that  the  derivaties  have  been  expressed  as  equivalent  distances  covered 
during  the  inter-sampling  period,  T,  and  that  the  predicted  acceleration  is 
unchanged  over  T,  since  a third  order  system  has  been  assumed.  Equation  (1)  can 
be  extended  or  contracted  for  higher  or  lower  order  tracking  filters. 

Matrix  notation  proves  very  useful  for  subsequent  analysis,  and  equation  (1) 
can  be  expressed  as: 


^t  +l/t  " \/t 


C2) 


where : 


A 

X 

1 

1 

1 

T.x 

; P = 

0 

1 

2 

O.ST^JC 

0 

0 

1 

(3) 


X is  a column  vector,  and  P is  the  prediction  matrix. 

The  Register  is  clocked  between  sampling  Instants,  to  hold  the  prediction  for 
the  next  correction  estimation  and  filtering  procedure,  at  time,  t+1. 


3.  THE  EXPONENTIAL  FILTER-PREDICTOR 


3.1  Derivation  of  filter  equations 

Consider  the  measurement  of  a physical  quantity,  x,  subject  to  noise. 

The  estimate  is  improved  if  the  average  of  several  measurements  is  taken, 
providing  the  noise  is  uncorrelated  between  seunples,  and  the  value  of  x does 
not  change  significantly,  during  the  sequence.  If  the  value  is  changing 
rather  slowly,  the  earlier  measurements  are  of  less  use  than  later  ones,  for 
estimating  the  current  value,  and  should  therefore  receive  less  weight  in  the 
average  calculation.  One  possible  estimator  gives  equal  weight  to  the  last 
N samples,  and  zero  to  earlier  ones.  This  requires  storage  of  the  last  N 
values,  so  that  subsequent  estimates  can  be  made.  The  exponential  estimator 
weights  all  samples  by  a factor  which  reduces  exponentially  with  the  time 
since  the  sample  was  taken.  Thus: 

K/t  “ “•  {*t  1 

The  factor,  a,  lies  between  one  and  zero,  and  determines  how  much  weight  is 
given  to  the  earlier  samples,  compared  with  later  ones.  The  coefficient, 
a,  outside  the  brackets,  normalises  the  estimate,  as  can  be  seen  by  putting 


Equation  (4)  can  be  expressed  recursively,  as  : 


t/t 


• ’‘t  * 

(5) 

\ 

t-l/t-l  * “•'•‘t  - 

(6) 

t/t-1  *■  “•  '"t  ■ ’’t/t-l' 

(7) 

Equation  (7)  follows  from  (6)  if  the  value  of  x is  assumed  to  be  unchanging. 
If  the  value  is  changing,  equation  (7)  intuitively  gives  a better  estimate 

than  (6)  (or  the  same  estimate,  if  derivatives  of  x are  not  used  in  the 

prediction  process).  Note  that  the  filtered  value,  is  the  predicted 

value,  plus  a times  the  difference  between  the  actual  measured 

value,  x^,  and  the  predicted  value.  Thus,  u times  the  measured  correction 

is  added  to  the  predicted  value.  ^ ^ 

By  an  extension  of  equation  (7),  the  filtered  values  of  T.x  and  O.ST^x 

are  obtained  by  adding  and  7 times  the  measured  correction,  x^^,  as 

corrections  to  the  predicted  values.  (This  is  the  reason  for  representing, 

X and  ft,  as  equivalent  distances  in  equation  (1).)  Thus  the  third  order 
exponential  tracking  filter  becomes: 
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t/t 


T.  X 


t/t 


O.Sfi  5, 


t/t 


*t/t-l  *ct 


*ct 


o.srt  ^ 


h/t-1  * *ct 


Or  in  matrix  notation: 


t/t 


^t/t-1  * *ct 


(8) 


(9] 


where : 


(10' 


3.2 


Filter-predictor  equations 

By  using  equation  (2)  to  eliminate  obtain  the  equation  defining 

the  operation  of  the  combined  filter-predictor,  via: 


\>l/t 


P.  X 


t/t-1 


+ P.F.  X 


ct 


(II 


This  can  be  rearrainged  to  give  a difference  equation: 


\+l/t  ■ ^t/t-1 


(P-I). 


't/t-1 


+ P.F.  X 


ct 


where,  by  evaluation: 


0 

1 

1 

a + 0 + y 

(P-I)  » 

0 

0 

2 

; P.F.  = 

0 * 27 

0 

0 

0 

y 

(ly 


(1 


3.3  Differential  approximation 

If  the  sampling  interval,  T,  is  sufficiently  small,  the  left  side  of 
equation  (12)  is  approximately,  T.  ^ . Taking  the  Laplace  transform  of 
(12),  and  using  this  approximation: 


T.  (s.X  - Xo)  * (P-I).  X ♦ P.F.  X 


where  Xo  represents  the  initial  conditions  of  the  vector  variable,  X. 
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By  rearraiiRcment , and  using  the  loft  Inverse  notion: 


i 


X (T.s  ♦ n.  1 - I' 


X + T.  Xo 

c 


By  evaluation: 


(14) 


f^TP 


I (T.s+l).r-P|“‘  = 


0 

0 


T.s 


2 + T.s 
2 T.s 
T".s’ 


(15) 


1 


|(T.s+l).I-Pr.  P-F  = T -s 


T*.  s^.  (a-t^  +7)  + T.  s . 03 +37)  +27 
T^.s^.  00+27)  + 2 T.s.7 
T^s^.7 


(16) 


I 


An  approximate  closed  relation  for  X can  be  obtained  from  the  inverse 
Laplace  transform  of  equation  (14) , providing  the  Laplace  transform  of 

X is  knovm. 
c 

3.4  Velocity  and  acceleration  lags 

Suppose  is  driven  (by  feedback)  to  cause  x to  follow  the  response: 


A 

X 


Xo  + Vo  .t  + 0.5-  ao  .t^ 


Then: 

X ^ Xo  ^ Vo  ^ ^ 

(17) 

s s^  s^ 

A .A 

By  inverting  the  approximate  transfer  function  relating  x^  and  x, 
from  equation  (14): 


X * 


T^s’  , 


X 

X 


f*.s’ 


(a-tfiyy)  + T.s.  (^+37)  + 27 


(18) 


The  steady  state  response  is  given  by: 


A 


A 


(♦OO 


Lim 

s^o 


(19) 
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From  equations  (17),  (18)  and  (19): 

Lim  j , T“.ao  , T.vo 

t>oo  ^ a 


(20) 


respectively,  for  the  third  order,  second  order  (7='o)  and  first  order  (7H3“o) 
filter-predictors,  showing  a finite  acceleration  lag  for  the  second  order, 
and  a finite  velocity  lag  for  the  first  order  case. 

A similar  conclusion  can  be  drawn  from  the  difference  equation  (12).  If 

A ^ • 

is  to  reach  the  zero  value,  the  change  in  X can  only  be  supplied  via  the 

matrix  (P-I).  Where  this  matrix  has  an  all  zero  row,  the  corresponding 

element  of  the  difference  vector,  A x,  must  also  be  zero.  Hence  we  require 
zero  changes  in  acceleration,  velocity  and  position,  respectively  for  third, 
second  and  first  order  filters,  if  no  lag  in  position  tracking  is  to  occur. 


4.  NOISE  FREE,  LINEAR,  TRACKING  LOOP  RESPONSE 


4.1  Exact  analysis 

A 

From  figure  1,  we  want  to  find  the  response,  ^ sequence  of 

inputs,  x^,  assuming  no  noise  and  a linear  discriminator  characteristic. 

Thus  5^  is  zero,  and  D is  a constant.  Examination  of  equations  (11)  and 

(13)  shows  that  the  coefficients,  a,  0 and  7 are  directly  multiplied  by  D, 
and  therefore  there  is  no  loss  of  generality  if  the  non-unity  value  of  D 
is  absorbed  into  them.  Therefore  D is  assumed  to  be  unity. 

Substituting: 


ct  t t/t-1 


in  equation  (12),  we  find  that: 


'^t+l/t  " \/t~l 


= (P-I  - P.F,).  X 


t/t-1 


+ P.F. 


(21) 


where : 


0 0 

.a~0.y 

1 1 

Fi  = 

0 0 

0 0 

; (P-I-P.F,)  » 

-0-27 

- 7 

0 2 

0 0 

(22) 


Again,  equation  (21)  can  be  evaluated  directly  by  a digital  computer. 

4.2  Differential  approximation 

By  analogy  with  the  procedure  of  section  3.3,  we  can  approximate  equation 
(21)  by: 


X { (T.s+l).I  - P P.F,j.-<.  { P.F.  X + T.Xoj 


(23) 
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where,  by  evaluation: 


I(T.s+l).I-P+P.F,j“‘  ' = 


|(T.s+l).I-P+P.F,j  .T'P’.F 


1 

E 


1 

S 


-T.s.  (P+27)-27 
-T.s.7 


T.s  2 + T.s  (24) 

t^.s*+T.s.  (a+/3+7)47  2 T.s+2o+^ 

-7  T^s^+T.s(a+<3+7)+^+27 


T^.s^.  (a+/}vy)  + T.s.  05+37)  +27 
T*. s’.  05+27)  + T.s. 27 
T’.s’.V 


(25) 


A = T^s’  + fls’ (a-H5+7)  + T.s.(p+37)  + 27  (26) 

Clearly,  A,  being  the  denominator  of  the  system  transfer  function,  is 
the  characteristic  function,  the  roots  of  which  determine  the  stability  of 
the  system. 

4.3  System  stability 

TTie  Routh  stability  criterion  (ref . 1)  for  a third  order  function,  reduces 
to  the  condition  that  the  product  of  the  two  inner  coefficients  exceeds  the 
product  of  the  two  outer  ones  (assiaming  all  are  positive).  Thus: 

(a+0+7)  (p+37)  > 27 


7 < 


P (a  * P) 

2 - 3a  - - 37 


(27) 


For  the  second  and  first  order  systems,  7=0,  and  the  systems  are  always 
stable. 

This  criterion  applies  only  to  the  differential  approximation.  Tbe 
actual  system  will  be  less  stable,  the  longer  the  period,  T,  due  to  phase 
shifts  produced  by  the  sampling  process.  If  T is  small  compared  with  the 
time  constants  associated  with  the  roots  of  (26)  (i.e.  a 1,  see  section 
5),  and  if  the  stability  margin  from  (27)  is  adequate,  further  analysis  is 
unwarranted.  Otherwise  an  analysis  based  on  the  z-transform(ref.2)  is 
required. 

4.4  Coefficients  for  critical  damping 

Being  of  third  order,  A has  three  roots,  one  being  real  and  two,  complex 
conjugates.  Classical  control  theory  suggests  that  optimum  performance  (in 
the  sense  of  best  compromise  between  noise  rejection  and  transient  decay 
time)  occurs  when  the  complex  roots  are  damped  somewhat  less  than  critically, 
but  that  performance  is  not  much  worse  at  critical  damping  (i.e.  when  the 
imaginary  part  of  the  conjugate  roots  is  zero).  For  simplicity  of  analysis 
we  will  examine  the  critical  case. 


mw- I K- 1 09«  (If  J 


If  the  time  constant  associated  with  the  third  root,  is  much  larger  or 
smaller  than  that  of  the  first  two  roots,  the  system  tends  to  approximate 
towarded  a first  or  second  order  system.  To  obtain  a typical  third  order 
system,  therefore,  we  choose  all  three  roots,  real  and  equal.  Thus: 


A = 

(T.s  ♦ P) 

a = 

1 - (1-P) 

0 

3P*  (2-p) 

2 

7 

Ph 

a+^+7  = 

3P 

0*37  = 

3p' 

(28) 


Naturally  (28)  satisfies  the  stability  criterion  (27)  with  a considerable 
safety  margin. 

4.5  Critical  damping  - second  order 

From  equation  (26),  with  7=0,  and  roots  real  and  equal 


A = (T.s  + pf  " 

a*0  = 2p 

(J  = p2 

a = p(2-p) 


(29) 


4.6  Optimum  damping  - second  order 

Benedict  and  Bordner(ref.3)  choose  as  the  value  giving  the  best 

compromise  between  noise  rejection  and  transient  decay.  This  gives  a 
1 

damping  of  ^ times  the  critical  value. 


A 


a 

+ — 
2 -a 


(1  + 


+ 


a 

2 -a 


(1 


P 


g 

2-g 


a 

0 


2^ 

1+P 

l^P 


(30) 
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4.7  Response  of  filtered  variable 

A 

Equation  (21)  gives  the  response  of  the  predicted  variable, 
Using  equation  (2),  the  response  of  the  filtered  variable  is: 


K/t  ■ " (P  - I " Fi  -P^^t-l/t-l  * ^-^t 


where,  by  evaulation: 


(P  - I - F,  .P) 


-a  1-a  1-a 


M -0  2-0 


-7  -7  -7 


As  before,  with  the  differential  approximation: 


X ^ {(T.s  + 1).  I - P + Fi  .Pj~l  jF.x  + T.xo  j 


where,  by  evaluation: 


{(T.s+l).I  - P + Fi  .Pj"*F 


j T*s^a  + T.s  03+7)  + 27 
T^.s^.P  + 2T.S.7 
T^s?7 


Naturally,  the  characteristic  function,  A,  is  unchanged.  Note  the 
trivial  difference  between  equations  (25)  and  (34),  with  7 ^0 

4.8  Velocity  and  acceleration  lags 

From  equation  (21): 

' Vi/t  - *t/t-i  "O’  - I - 

- P.F.  Xj  - (P  - I - P.P.). 


where  X is  the  column  vector  representing  the  input  variable,  and  X^^j  » 
P^.X^  (analogous  to  equations  (3)  and  (2)).  Assuming  steady  state 

conditions  are  reached  after  a sufficiently  large  time,  the  left  hand  side 
becomes  zero.  Thus: 


(X-X) 


t/t-1 


- p + p.pir*  [(Px-p)-\/t-i  * *Vt-i  ■ ’’•v] 


(I  - P + P.F,)-*.  (P^  - P).x 


t/t-1 


since  the  second  term  within  the  compound  brackets  is  zero. 
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For  a second  order  processor  with  constant  acceleration  input: 


1 1 0 


P = 0 1 0 ; 7 = 0; 


0 0 0 


1 1 1 


0 12 


0 0 1 


By  evaluation: 


0 0 2/0 


Cl  - P + P.Fir'.  (Pj^-P) 


0 0 


0 0 1 


Whence : 


2 . . TO  , 


x-x  = 0-.  O.SaoT^  = ^ 


• ^ 2a*0  O.S.ao.t^  _ ao.T  f2a-W3: 

x-x  ~ ^ • T ~ 20 


For  a first  order  processor  with  constant  velocity  input: 


By  evaluation: 


1 0 


0 0 


; 0 = 7 = 0;  P, 


Cl  - P + P.F,n  XP^  - P) 


1 1 


0 1 


0 i 

a 


0 1 


Whence : 


1 

x-x  = ;r.T.Vo 


5.  TYPICAL  TRANSIENT  RESPONSES 


Figures  2 to  6 inclusive  illustrate  the  transient  decay  of  the  position, 

A 

velocity  and  acceleration  corrections  C*~X) » following  step  changes  in  position, 
velocity  and  acceleration  of  the  measured  variable.  Equation  C21)  is  evaluated 
directly  by  digital  computer,  for  values  of  a,  0 and  7,  from  Table  1,  being 
appropriate  to  equal  root  critical  damping  for  first,  second  and  third  order 
trackers,  and  to  optimum  damping  for  a second  order  tracker.  These  are 
labelled,  1,  2C,  3,  2q  respectively. 
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The  curves  are  normalised.  The  abscissa  is  calibrated  in  terms  of  , 

T 

i.e,  time  normalised  to  the  basic  time  constant  of  the  filter  (T  = — ) • The 

T " 

normalised  sampling  interval  Cp  = P)  is  0.02  for  the  continuous  curves,  which 

are  indistinguishable  from  curves  for  the  limiting  case  of  p = 0.  The  crosses 
(shown  only  for  the  first  and  third  order  trackers)  mark  the  responses  for  P = 0.2. 

The  locus  of  the  response  for  p = 0.2  is  reasonably  well  approximated  by  the 
limiting  curve  (P  = 0),  differing  mainly  in  greater  overshoots  and  slightly 
shorter  decay  times.  The  first  order  tracker  has  a good  recovery  from  a position 
step  (figure  2);  tolerates  a velocity  step  with  a lag,  providing  the  lag  error 
does  not  exceed  the  drop  out  range  of  the  discriminator  (figure  3) ; and  rapidly 
loses  lock  under  a prolonged  acceleration  step  (figure  4).  It  gives  the  best 
noise  rejection  on  the  position  estimate  (see  Section  6,  and  Table  2). 

The  second  order  tracker  handles  a velocity  step  without  lag,  (figure  3)  and 
an  acceleration  step  with  lag  (figure  4).  In  all  aspects  the  optimum  damping 
case  is  superior  to  the  critically  damped  case  (i.e.  reduced  lag  and  recovery 
time).  It  is  comparable  with  the  first  order  system,  under  a step  in  position. 

It  gives  a velocity  estimate,  but  with  a noisier  position  estimate  (Table  2). 

The  third  order  system  is  inferior  in  all  respects,  except  for  providing  an 
acceleration  estimate  and  no  acceleration  lag.  The  recovery  time  could  perhaps 
be  improved  by  separating  the  roots  of  the  characteristic  equation,  to  reduce 
the  long  overshoot  (figure  2). 

The  lags  check  with  the  results  calculated  in  Section  4.6,  and  listed  in 
Table  1.  For  a given  time  constant,  t,  and  damping  factor,  the  lags  are  only 
slightly  dependent  on  the  normalised  sampling  time  (p) . 

6.  NOISE  REJECTION  BY  LINEAR  TRACKER 

6,1  General  theory 

We  consider  the  linear  exponential  tracking  loop,  illustrated  in  figure  1, 
with  D = 1,  as  discussed  in  Section  4.1,  In  this  case,  the  external  signal, 
is  arbitrary,  and  will  be  assumed  to  be  zero.  The  noise,  6^,  is  non-zero. 

(Note  that  it  does  not  include  output  noise  fed  back  via  since  this 

latter  terra  appears  explicitly  at  Thus  8^  can  be  considered  to  act 

at  X . Hence,  we  wish  to  obtain  the  standard  devi#«’'isri  of  the  predicted 

^ A 

variable,  ^s  a function  of  the  standard  deviation  of  the  input 

variable,  x^,  after  steady  state  conditions  have  been  achieved. 

From  equation  (21): 

Vi/t  " Vt-1  * \ 

Writing  the  first  and  second  rows,  explicity: 


"t+l/t 


= qn.  ♦ 


+ r,  .x. 


‘t+l/t  ""  *t/t-l  * ‘^”*^-\/t-l  * 


♦ rj  . X. 


Therefore : 


Lim  \ jq- 


N*’‘» 


Lim 

N^oo 


t=0 

N 


N 1 ^ 


t=0 


1 A 1 ^ 


t=0 


t=0 


It  is  clear  that  is  uncorrelated  with  all  other  terms  on  the  right  hand 

side,  since  they  are  based  on  data  at  (t-l).  Under  steady  state  conditions, 
the  limit  on  the  left  hand  side  has  the  same  value  as  the  corresponding  limit 
product  term  on  the  right  hand  side  (ignoring  the  q coefficients).  Thus,  for 

M 

a system  of  order  M,  we  have  (M  + C2  •)  equations,  expressable  in  the  form: 


V. 


Var  (x) 

'n^ 

Var(T.x) 

• 

V 

Covar(x,T. x) 

a 

ri  . T2 

ri  .T3 

.Var (x) 


(41) 


where: 


(J-qii)  -q?j  -qiJ.  . . 
-qJi  (l-qJi)-qjj.  . 


•qj^  j *2qii.qii  -2qii.qi], 
*q*2M  I ■2‘\n-qii  *2qii.qj3. 


-H, 


• • • • 

• • • • 


Ml 


-qii-qji  -qij-qii  .•••  -q,u-qjul(*'‘ii' '(qi i -qjs^qi  j .qn ) . 


'1M'^2M 

■ -‘’iM-^SM 


■^‘*1(M-1)’’2M* 


r*’(M-l)''’Mr  • • • '''(M-1)M'%wI'^'’(M-1)1‘*'m2*‘'(M-1)2-%«‘ 


(42) 


- 13  - 


»fRE-TR-1854(W) 


Thus  the  steady  state  variances  and  covariances  of  the  elements  of  X,  are 
obtained  by  multiplying  equation  (33)  by  the  left  inverse  of  V,  a square  matrix 
M 

of  order  (M  ♦ Cj  ).  Evaluating  the  inverse,  algebraically,  is  hardly  feasible 
for  a system  of  third  order,  or  higher. 

6.2  Prediction  covariance  - second  order 

For  a second  order  system,  M = 2,  and  the  matrix,  V,  is  of  order  3.  From 
equations  (21),  (22)  and  (13): 


Q = P - P.F, 


1 - CL  - 0 


* P.F  = 


a * P 


By  evaluation: 


(2  - a -^)(a  + (?)  -1  -2(1  - a . p) 


00  (4-2a^5) 


P(1  - a - p) 


-1  a + 2P 


0^  (2-a) 


And  by  evaluation  of  equation  (33): 


-a^  ~2CL^P  -oP^ 
*200  *2a^  *20 

CL  (2-n)  *P  (1-a) 


-20^  (l-a) 
-200  (2-a) 


VarCS)  = ° * . Var(x) 

2 - a .012 


Var(T.x)  = ^ (-  — . Var(x)  (44 

2 - a .012 

Covar(x,  T.l)  = a 

6.3  Filtered  covariance  - second  order 

For  the  filtered,  rather  than  predicted  variable,  we  use  equations  (31), 
(32)  and  (10): 


I 

I 


imE-TR-1854  00 


> (P  - I t .1*) 


R = F * 


1-a  l-a 


l-fi 


a(2-a)  -(l-a)’ 


-2(l-a) 

20(1-0) 


0(1-0)  .(i.a)(l-0)  0+20-200 


<10(4-20-0) 


0(2-00)  2-5a+4a’-a^ 


20(1-0)2 


02(2-0)  a2  (2-a)+20(l-a)  20(0-2a+a^)' 
-0^  (l-a)  2o-3o2 +d2+«0-0  0(4o-2o2-0) 


Var(x)  = 


VaT(T.x)  = 


o + 0/a  - 30/2 
2 - o - 0/2 


. Var(x) 


i 0/2  • 


Covar(x,  T.x)  = • Var(x)  (48) 

Note  that  these  values  are  marginally  smaller  than  those  for  the  predicted 
variable  (except  for  identical  velocity  variance)  as  would  be  expected  from 
general  statistical  considerations. 

6.4  Covariance  - third  order 

Haber (ref. 4)  has  succeeded  in  deriving  an  exact  relation  for  the  covariances 
of  the  third  order  filtered  variable,  for  the  critically  damped  case,  (although 
there  is  an  error  in  one  of  them).  Correct  results  are  listed  below.  He 
uses  an  ingenious  matrix  factorization,  which  reduces  the  algebraic  manipul- 
ations to  a reasonable  level.  (Note  that  Haber's  "o"  is  the  "p"  6f 
this  report.  Section  4.4). 


Var(x)  = 


pf  66  - 186P  + 202p2  - 100p2  + I9p*  j.  Var(x) 

(2-p)* 


Var(T  X)  - ^2  - 148p  49p2i Var(x) 

2(2-p)* 
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Var(0.5T^.  x)  = 3p-^  (51) 

2(2-p)® 

6.5  Summary  of  noise  rejection 

A 

Table  2 gives  approximate  values  of  the  ratio  of  variance  of  X to  the 
variance  of  x,  derived  from  equations  (46)  to  (51) , in  terms  of  the  basic 
real  root,  p,  for  critical  and  optimum  damping,  as  defined  in  Sections  4.4, 

4.5  and  4.6. 

As  would  be  expected  from  the  statistical  theory  of  variance,  the  trackers 
with  more  degrees  of  freedom  (higher  order)  show  the  smallest  reduction  in 
noise.  Thus  the  system  order  should  be  as  small  as  consistent  with  accept- 
able lags  with  likely  targets. 

Note  that  in  each  case,  the  effective  time  constant  r of  the  tracking  filter 
is  given  by  T/p. 

7.  EQUIVALENT  CONTINUOUS  DATA  TRACKER 

Figure  7 shows  a typical  third  order  continuous-data  tracking  filter-predictor 
consisting  of  a velodyne  integrator  (loop  A3  - B)  preceded  by  two  integrators,  Aj 
and  Ai ,.  stabilised  with  resistors  in  series  with  the  capacitors.  These  provide 
the  well  known  proportional  plus  integral,  error  feedback.  In  the  ideal  case, 
with  infinite  amplifier  gains: 


- xo  . - 


x .dt  + X 
c c 


= ^ + ’’2  ^.dt  + X 


Xo  + — x.dt 

" 3 


where : 


1 2 = R . C2 


Taking  the  Luplace  transform  of  each  equation; 


— S'  — 1 

A Xo  * n ^ ^ 

* = T • ■‘c-o  * rf;’ 


r '5- 


WRE-TR-18S4(W) 


16  - 


Multiplying  the  first  and  second  equations  by  O.ST*  and  T,  respectively,  to 
preserve  the  same  form  of  vector  used  previously  (T  has  no  other  significance, 
and  cancels  in  practice  when  equations  are  solved): 


1 

1 

n 

i\ 

1 

T.s.Tj 

U 

A 

u 

0 

1 

l+S.Tj 

0.5T.STJ 

T. 

.X 

= 'b. 

s 

+ 

0 

0 

0 

1 

0.5T*  (l+s.Ti  ) 

S.T, 

.X 

c 


By  multiplying  through,  by  the  left  inverse  of  the  left -most  matrix: 


1 

1 

l+S.Tj 

1 

(l+S.T,)(l+s.r2) 

T.  s .r3 

0.5T*.S*.Ti  .Tj 

S^T, .TjTs 

K 

X = 

0 

1 

1+S.T, 

0.5T.s.ri 

A 

Xo 

. — + 
s 

(l+S.T,)(l+S.T2).T.S.r3 

0 

0 

1 

(1+s.r, )0.5T^s?r2 .r3 

(53) 

A comparison  of  equations  (47)  and  (14),  (15)  and  (16)  shows  the  basic 
similarity  of  form,  although  minor  differences  exist.  The  transfer  function 

between  x and  can  be  made  identical  by  making: 


T 

a+/J+7 

’’i  .’’i 

^ P*3y  t y/0-yf  - 80.7 

47 

P*3y 

r,+T, 

a T.  

27 

.Tj 

, ^2  a^*y 

27 

.fj  .T3 

= T*/27 

For  a second  order  filter. 

7 = 0 and  Ti  » “>  , whence: 

^2 


T. 


(54) 


(55) 


T,4T, 


KK 
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Table  3 shows  the  values  of  , Tj , Ts  for  the  first  second  and  third  order 

filters  considered  previously.  It  happens  that  the  configuration  of  Figure  7 

is  not  suitable  for  producing  a third  order  filter  with  three  equal  real  roots, 
as  the  required  values  of  r,  and  r,  are  complex.  The  values  given  in  the  table 
gives  roots  P,  P and  0.25P.  Thus  one  root  gives  four  times  the  time  constant  of 
the  other  two. 

The  equations  (54)  combined  with  equations  (46)  to  (48)  provide  a means  of 
estimating  the  noise  rejection  of  the  continuous  data  filter,  fed  with  data  with 
correlation  time  less  than  T.  Thus,  for  a pulse  radar,  the  time,  T,  corresponds 
to  the  inter-pulse  period. 


8.  CONCLUSION 

Basic  equations  have  been  derived  for  the  third  order  exponential  tracking 
filter-predictor  operation  and  noise  rejection.  Normalised  curves  are  given 
for  the  transient  decay  following  position,  velocity  and  acceleration  steps, 
for  critically  damped  first,  second  and  third  order  filters,  and  for  an  optimally 
damped  second  order  filter.  These  apply  exactly,  when  the  sampling  interval 
is  small  compared  with  the  filter  time  constant,  and  are  fair  approximations  for 
ratios  of  up  to  0.2.  Performance  for  other  values  can  be  obtained  exactly  by 
digital  computer. 

The  second  order  optimally  damped  filter  gives  best  results  providing  the 
acceleration  lags  can  be  tolerated.  This  depends  on  the  expected  variation 
with  time  of  the  primary  variable. 

A third  order,  continuous  data  (analogue  type)  of  tracking  filter  has  also 
been  examined,  and  equivalence  relations  with  the  exponential  filter  have  been 
derived.  This  permits  direct  calculation  of  the  noise  rejection  for  analogue 
trackers,  based  on  the  decorrelation  time  of  the  input  correction  data. 

In  both  continuous  and  sampled  digital  systems,  the  noise  improvement  is 
proportional  to  the  equivalent  number  of  uncorrelated  correction  estimates 
received  per  filter  time  constant;  the  lower  the  order  of  the  filter  the  better 
the  improvement. 
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NOTATION 

discriminator  function  (Figure  1) 
column  vector  (equation  10) 
column  vector  (equation  22) 
unit  matrix 
order  of  filter 

prediction  matrix  of  filter  (equation  3) 
prediction  matrix  of  input  variable 

general  iterative  matrix 
general  iterative  column  vector 
intersample  period 

covariance  generating  matrix  (equation  42) 
input  variable  vector 
initial  condition  vector 
correction  vector 

initial  acceleration  of  input  variable 

imaginary  operator 

element  of  Q 

element  of  R 

Laplace  operator 

time 

initial  velocity  of  input  variable 
element  of  X 

initial  position  of  input  variable 
correction  in  position  (element  of  X^) 

characteristic  equation  (26) 

filter  parameter 

filter  parameter 

filter  parameter 

noise  sample 

root  of  A 

filter  time  constant  (T/p) 

time  constants  of  analogue  tracker 

differentiation  with  respect  to  time 

estimated  value  of  

Laplace  transform  of  

suffix  indicating  value  of  t,  predicted  at  time  (t-1) 
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TABLE  1.  ESTIMATE  STEADY  LAGS  DUE  TO  TARGET  MOTION 


TABLE  2.  NOISE  REJECTION  PERFORMANCE  OF  TRACKERS 


Var(T.x)/  . Var  (0 . 5T».  x)  / 

Var(x) 
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TABLE  3.  PARAMETERS  OP  CONTINUOUS  FILTER- PREDICTOR 
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Figure  5. 


Velocity  response  to  position,  velocity  and  acceleration 
steps 
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Figures  6 and  7 


Figure  6.  Acceleration  response  to  position,  velocity  and  acceleration 
steps 


Figure  7.  Typical  analogue  tracking  filter -predictor 
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